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Neutral hydrogen is ubiquitous, absorbing and emitting 21 cm radiation throughout much of 
the Universe's history. Active sources of perturbations, such as cosmic strings, would generate 
simultaneous perturbations in the distribution of neutral hydrogen and in the Cosmic Microwave 
Background (CMB) radiation from recombination. Moving strings would create wakes leading to 
21 cm brightness fluctuations, while also perturbing CMB light via the Gott-Kaiser-Stebbins effect. 
This would lead to spatial correlations between the 21 cm and CMB anisotropies. Passive sources, 
like inflationary perturbations, predict no cross correlations prior to the onset of reionization. Thus, 
observation of any cross correlation between CMB and 21 cm radiation from dark ages would 
constitute evidence for new physics. We calculate the cosmic string induced correlations between 
CMB and 21 cm and evaluate their observability. 



I. INTRODUCTION 

Neutral hydrogen emits and absorbs radiation at 21 cm 
in its rest frame throughout the history of the Universe. 
The brightness of the 21 cm line observed today from 
various redshifts is determined by the spatial distribu- 
tion and peculiar velocity of the neutral hydrogen that 
emitted it. The radiation from any particular redshift has 
small angular variations, similar to those in the Cosmic 
Microwave Background (CMB) radiation. There is also 
variation along the line of sight, with a correlation depth 
of approximately 10 Mpc PQI2]- I n the absence of physics 
that violates free streaming, one does not expect a corre- 
lation among slices in the line-of-sight direction that are 
separated by more than this correlation length. Since 
observations can measure the redshift of 21 cm emission 
precisely, we should be able to extract information about 
this line of sight variation from 21 cm observations. 

It is this third dimension of information encoded in 
21 cm radiation that can make it a powerful new tool to 
learn about fundamental physics. One path forward is to 
look for physics that violates free streaming and gener- 
ates correlations among distinct redshift shells [2i . Sev- 
eral sources of such a signal have been identified, though 
they are limited to the epoch during and after reioniza- 
tion [SHU- Though passive perturbations, such as those 
seeded by inflation, cannot generate such correlations, 
active sources of perturbation could. Active sources of 
perturbation include cosmic strings and other topological 
defects, such as textures [7], as well as the field gradients 
that are left behind after global phase transitions [8] . Of 
these, cosmic strings are the most thoroughly studied and 
have perhaps the richest phenomenology [9j , so they will 
be our focus. 

Moving strings generate density wakes in the material 
through which they pass [TO] , including neutral hydrogen. 
As a result, the brightness of the 21 cm emission, which 
depends on the density and the velocity of the emitters, 
will be directly perturbed in the string's wake [IT. The 
same moving strings also perturb the background CMB 



light, monotonically shifting the spectra and creating line 
discontinuities in the temperature screen of the sky — a 
phenomenon known as the Gott-Kaiser-Stebbins (GKS) 
[l2j H3] effect. GKS is a special case of the Integrated 
Sachs- Wolfe (ISW) effect [H] caused by time-varying 
gravitational potentials sourced by moving strings. Thus, 
there will be some correlation between string-sourced 
21 cm brightness fluctuations at any redshift and a part 
of the CMB temperature anisotropy. Strings will also 
induce some cross-correlation in 21 cm fluctuations com- 
ing from different redshifts. However, as we discuss in 
Section [nj we expect this effect to be relatively small. 

Precision measurements of the CMB temperature 
anisotropy limit the contribution of cosmic string pertur- 
bations to be less than about 10% of the total over the 
range of scales covered by WMAP [T5] . As such, the bulk 
of the Universe's anisotropy comes from the primordial 
inflationary fluctuations, for which the cross-correlation 
between sufficiently-separated shells is zero. The fact 
that any string contribution to primordial anisotropy 
must be small means that the intrinsic signal for which we 
are searching will be faint. However, since there is noth- 
ing in the standard ACDM model that could generate 
such cross correlations before the epoch of reionization, 
we can be sure that discovery of such a cross correlation 
would be a sign of new physics. 

To examine the prospects for detecting these string- 
induced correlations, we will evaluate the signal-to-noise 
density of this effect and the volume over which the sig- 
nal can be found. Our aim is to identify the optimal way 
of looking for stringy phenomena in the 21 cm radiation 
from the dark ages. This effect is present over a vast 
volume — from z < 100 until the onset of reionization — 
which may make a detection of the CMB/21cm correla- 
tion feasible under optimistic observing scenarios. How- 
ever, because of uncertainties in reionization physics, we 
exclude a large region (z < 20) from our present study. 
We will also leave the calculation of 21cm x 21cm cross- 
correlation from different redshifts for future work. 

In this paper we adopt the string network model of 
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[TBI IT?] , which was previously implemented in CMB- 
FAST pi] and made available publicly as CMBACT [TJ]. 
For this work, we have implemented this string model in 
CAMB [5D]. Our code will be made publicly available 
soon [2"T1. 



II. 21 CM/CMB CORRELATION INDUCED BY 
COSMIC STRINGS 



Cosmic strings are relativistically moving linear topo- 
logical defects that are remnants of the high energy early 
Universe. As they move, strings affect matter by kicking 
it into over-dense wakes and perturb light by shifting its 
spectrum via the GKS effect. Strings form in the early 
Universe and remain active through all epochs, spreading 
their effects throughout the Universe's history. 

The effect of cosmic string networks on the CMB power 
spectrum has been well studied. Strings affect the CMB 
in both of the ways described above: their wakes gener- 
ate a single peak at an angular scale determined by the 
strings' correlation length (< c/H), and strings in the 
foreground of the CMB generate an interlacing network 
of linear discontinuities through the GKS effect. The fact 
that these signatures have not been observed provide a 
bound on the density and tension of cosmic strings [15] . 

These effects are also present during the dark ages, 
when the neutral hydrogen is releasing 21 cm radia- 
tion. Moving strings produce wakes and thus seed hydro- 
gen density contrasts that grow and generate observable 
variations in the brightness temperature. Interestingly, 
though, the GKS effect is much less observable for 21 cm 
radiation. Strings shift the spectrum of light that passes 
them, so their effect on passing light is only observable 
when there is a noticeable change in the known back- 
ground source, such as a shift in a spectral peak or dis- 
continuity in a smoothly varying brightness. However, 
in the presence of a flat spectrum of background light, 
the string is invisible, since a lateral shift of a horizon- 
tal spectrum results in no detectable change. In the dark 
ages, strings will be swimming in a background of contin- 
uously renewed radio waves, with a very weak gradient in 
brightness and spectrum. Thus, the strings' GKS effect 
on 21 cm radiation will be considerably less observable 
than the same effect applied to CMB photons 1 . 

To calculate these string-sourced effects quantitatively, 
we must give a formal characterization of the string 
perturbations. The most important difference between 
string-sourced perturbations and those from inflation, be- 
sides the active/passive distinction, is that strings source 
vector mode perturbations. For strings, vector modes are 
as large as scalar modes, and they also generate much of 
the novel phenomenology of string networks [22j. Vec- 
tor modes are usually neglected since they decay in sys- 
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terns described by General Relativity and ordinary mat- 
ter. Strings, however, actively source vector perturba- 
tions faster than they can decay away. We will reproduce 
the most relevant equations for vector mode perturba- 
tions in the Appendix. 

Calculating CMB and 21 cm observables also requires 
a string model and a code for tracing the history of light 
streaming to us now from these early times. To this end, 
we ha ve in corporated a model for strings [111 [T71 [T5] 
(see § II A for a description) into the 21 cm extension 
of CAMB [20 ; this code will be publicly released (2TJ. 
Our major addition to the underlying CAMB code is the 
inclusion of vector mode fluctuations. Details about how 
string-sourced vector modes enter the Boltzmann equa- 
tion can be found in Refs. [TBI H3 EQl [23] , and will be 
more comprehensively explained in |21j . 

Now we can ask what physics will control the shape 
and amplitude of string-sourced two-point functions in 
the 21 cm fluctuations and in the cross correlation. The 
two physical effects controlling the shape of these spec- 
tra are 1) the intrinsic brightness of the 21 cm photons 
and 2) the correlation length of the string network. The 
brightness temperature sets the amplitude of the spec- 
trum. It, in turn, is determined by the difference between 
the 21 cm spin temperature and the CMB temperature 
at that epoch, which are mismatched over a long period 
during the dark ages [2]. It is worth noting that the 
brightness temperature is negative throughout most of 
the dark ages. That is, the neutral hydrogen is absorbing, 
not emitting, radiation. This shows up as an observable 
because we can infer what the flux of long wavelength 
CMB photons should be, and hence detect the dimming 
due to neutral hydrogen. Since we are looking for small 
variations in the brightness temperature, a large intrin- 
sic brightness makes detections easier. The anisotropy 
of 21 cm brightness is chiefly caused by variations in the 
density and the peculiar velocity of the neutral hydro- 
gen, both of which are sourced by cosmic strings. These 
effects and others are more comprehensively discussed in 
Ref. [2D]. 

The correlation length of the string network sets the lo- 
cation of the peak in £-space, the angular scale. Since this 
correlation length is some fraction of the Hubble scale, 
the spectrum from dark ages strings peaks for multipolcs 
I < 100. This peak in the two-point function corresponds 
to the size of the wakes generated by the string network 
at that epoch. Hence, the peak signal migrates to lower-f 
as the redshift decreases since the network is scaling with 
the horizon, and the horizon is growing. This trend can 
be seen in our Fig. [TJ a plot of the CMB-21 cm cross- 
correlation Cf z (z,£) from a string network, which we 
have generated with our code. 



A. The string model 

The string code we used is based on the segment model 
[161 117j . In this model, the string network is represented 
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FIG. 1: CMB-21 cm cross-correlation power spectra 
(Cj z [/iK] 2 ) generated by a network of strings with tension 
Gfi — 4 x 10 -7 . The amplitude of the signal scales with the 
brightness temperature, peaking at a redshift z ~ 55. The 
general shape — anti-correlation with a single broad peak — is 
similar for all redshifts, though the £ of maximum correlation 
moves to larger scales at later times / lower redshifts. This 
pattern comes from the string network's scaling behaviour: 
strings seed fluctuations at a fixed fraction of the horizon size 
at the epoch in which the seed is created. This corresponds 
to larger scales (smaller i"s) for later times, i.e., for lower 
redshifts. A contour at —4 x W~ 6 [/iK] 2 is shown. 



by a collection of uncorrelated straight string segments 
moving with uncorrelated, random velocities. There are 
two fundamental length scales in such a model: £, the 
length of a string segment, which represents the typ- 
ical length of roughly straight segments in a full net- 
work; and £, the typical length between two string seg- 
ments, which sets the number density of strings in a given 
volume (N s oc l/^ 2 )- The model also tracks the root- 
mean-square (rms) velocity of the string segments. These 
parameters are deduced by assuming network scaling — 
i.e., the string energy density tracks the dominant back- 
ground energy density in the radiation and matter eras. 
Numerical simulations of networks have confirmed that 
these assumptions accurately describe network physics. 
The positions and orientations of the segments are drawn 
from uniform distributions in space and on a two-sphere, 
respectively. The model's parameters have been cali- 
brated to produce source correlation functions in agree- 
ment with those in Ref. [24]. Because of the approxima- 
tions built into the string segments, the model can only 
produce two-point correlation functions. It cannot be 
used for making simulated string maps or other more so- 
phisticated observables; though it could, with some mod- 
ification, generate three-point functions |25j . 

On the cosmological scales we consider, finer details of 
the string evolution do not play a major role. It is the 
large-scale properties — such as the scaling distance, the 
equation of state (wiggliness), and the rms velocity — that 
determine the shape of the string-induced spectra. Fur- 
ther details of the computational methods necessary to 
compute the 21 cm signal will be included in a subsequent 



publication [2Tj . 

B. Cross correlations 

Let 5T z (x,n) denote the anisotropy in the 21 cm 
brightness temperature emitted at redshift z and ob- 
served at a spatial position x coming from the propa- 
gation direction n. Similarly, the CMB brightness tem- 
perature is characterized by 5T CMB (x', n'). For com- 
putational purposes these quantities are expanded into 
normal modes |23) 

5T X (x, n) = [ ^ Y, E 6T £nik) GT , (1) 

J ^ ' t rn=-2 

where k is the Fourier conjugate to x, X S {z, CMB} 
and STf m are the associated multipole moments in the 
basis 

° T = H) ' \I2IT1 Yr (fl) 6XP O*'*)' (2) 

In our analysis we use the angular correlation function 
between quantities X and Y as measured by an observer 
at x — 0, defined as 

C XY {Q) = ( < 5T z (n) ( 5T CMS (n')>fr,v=cose 

1 OO 

= r J2( 2£+1 ) c * Yp e( cos@ ) ( 3 ) 

where Pi are the Legendre polynomials, and Cf Y is the 
anisotropy spectrum calculated from the associated mul- 
tiples 8T x {k) PH [2D] as 

Cf Y = (4-) 2 / d*k J2 5T? m 5T? m . (4) 

In the above notation, m = corresponds to scalar 
modes, m — ±1 labels vector modes of right- and left- 
handed vorticity, and m = ±2 denotes the tensor modes. 
See the Appendix for an expanded discussion of our cal- 
culational method. 

III. SOURCES OF NOISE 

Having formed the cross correlation spectrum 
(Eq. Q), a density distribution of signal-to-noise (S/N) 
can be formed in redshift and I space as 

Az is the width of the redshift shell set by the chosen 
frequency interval and / s k y — 0.8 is the fraction of ob- 
servable sky. The total signal-to-noise (S/N) is then 
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where the sum is over the harmonic multipoles £ and the 
(independent) observed redshift slices z% of width Azi, de- 
termined by the observational setup. Each multipole can 
receive contributions from the inflationary background 
Cf d , the active sources Cf tr , and astrophysical and tele- 
scope noise , so we can write 

Cf Y [total] = Cf XY + Cf XY + C XXY . (7) 

Strings are the only source of the cross-correlation, so 
the numerator in Eq. ^ receives only one contribution. 
The denominator, however, receives a large contribution 
from inflationary perturbations that are at least ten times 
the size of the string power spectrum. There is also a 
large contribution from the detector noise which is dom- 
inated by the galactic emission at the low frequencies we 
need to probe the high redshift 21 cm brightness tem- 
perature. The detector noise is uncorrelated between the 
CMB and 21 cm frequencies and is given by 



where T sys is dominated by the sky temperature and, for 
regions away from the galactic plane, can be approxi- 
mated as [2] 

T svs (v) = 150K ( ) " . (9) 

y V ' \ 180MHz J y ' 

In the previous expressions Av is the width of a frequency 
bin, which is set by the experiment, i ODS is the observing 
time, / covor = iVdishAiishMtotai is the fraction of the total 
area enclosed by the experiment and characterizes the 
total collecting area, and i m ax{v) = 2irD a _ rray v / 'c, where 
-Darray is the diameter of the array and c is the speed of 
light.' 

Foregrounds include galactic and extragalactic point 
sources, as well as the galactic and extragalactic free- 
free emission, all of which are expected to induce correla- 
tions between different frequencies [2H HZ] . In this study, 
however, astrophysical sources will have a negligible ef- 
fect since we are only including 21 cm cross-correlations 
with CMB photons, and these signals quickly decorrelate 
with increasing frequency separation [2|. Furthermore, 
these signals dominate on very small scales while string- 
induced correlations reside on large scales, t < 100. 

We fold the observational parameters into a single ob- 
servational ambition parameter, x, defined as 

x = t ohs f covcr D 2 , (10) 

with units of [x] = [years] x [km] . There is always a 
practical trade-off between making f COVCI large and mak- 
ing D large. For our signal, though, it is clear what we 
should choose. The main reason for going to large D 
is to resolve very small angular scales. However, string 
sourced cross correlations reside only on large angular 
scales. Hence, we can make due with a relatively small 
baseline D, which should let us make / COVO r nearly unity; 
i.e., a "moderately" sized telescope or array with nearly 
a complete covering area would be ideal for searching for 
strings. 



IV. RESULTS 

We have calculated cross correlation power spectra 
from redshifts throughout the dark ages. In producing 
our results, we have fixed the observing bandwidth to 
Av = 0.5MHz; however, varying this does not change our 
result. Although increasing the bandwidth diminishes 
the autocorrelator signal (noise "iV"oc l/Av) and, thus, 
increases the signal-to-noise in each redshift bin, one ob- 
serves fewer independent redshift slices as a result (signal 
"Stot'oc 1/Az/) so the overall signal-to-noise remains un- 
changed to first order. We have numerically verified this 
for bandwidths in the range 0.01MHz < Av < 2MHz. 

We have excluded the epoch of reionization (z < 20) 
from our signal-to-noise results. There are several rea- 
sons for this. The main reason is that reionization physics 
are still too uncertain for reliable calculations to be done, 
and the details of reionization, once begun, will dom- 
inate the 21 cm physics. Moreover, if cosmic strings 
are present during reionization, they will likely play a 
role in speeding up the process of reionization through 
their early generation of non-linear overdensities [28l [29] . 
Reionization itself will also create cross-correlations be- 
tween the CMB and 21 cm radiation [3HS]- Though it 
is possible that the string and reionization signals could 
be disentangled if they are both present, we restrict our- 
selves in the present study to redshifts where strings are 
the only possible source of cross-correlation. 

At redshifts z < 7 reionization is complete and the sur- 
viving neutral hydrogen only exist in islands that are pre- 
sumably correlated with large scale structure. In princi- 
ple, we can continue calculating the string sourced cross 
correlation between the matter distribution (traced via 
21 cm or any other observable light) and the CMB. How- 
ever, there will be competing sources of cross correlation. 
For instance, the re-scattering of CMB photons off ion- 
ized particles will induce a correlation of CMB with the 
large-scale plasma distribution |30j . Eventually, at lower 
redshifts, there will be a correlation between large scale 
structure and the ISW distortion of the CMB caused time 
varying potentials when Dark Energy becomes dynami- 
cally important. Nonetheless, we have made an approxi- 
mate calculation of what a hypothetical radio survey at 
these redshifts might see in the presence of strings, ex- 
tending an analytical distribution of large scale structure 
tracers used in [3T] to the range < z < 7. The resulting 
theoretical signal-to-noise for the cross correlation was 
too small. This is not surprising since string wakes that 
are sourced at such late times have not had time to grow. 
By contrast, the potential wells that attract the ionized 
particles which then re-scatter CMB have been growing 
through gravitational instability for some time since en- 
tering the horizon. Hence, we do not include the z < 7 
range into our signal-to-noise results. 

Let us turn to our results, shown in Fig. [2] which begin 
with analysis of a hypothetical Universe where inflation- 
ary perturbations are not present and cosmic strings are 
the only source of both CMB and 21 cm anisotropy. We 
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FIG. 2: Distribution of signal-to-noise squared (S/N) 2 from cross-correlation studies between CMB and 21 cm photons in the 
presence of a string network. The observing bandwidth has been fixed to Av — 0.5MHz, which determines the redshift bin 
width | Az\ = (1 + z) 2 Av/vq, vq being the rest frequency of the 21 cm radiation. Signal along horizontal slices is modulated by 
the brightness temperature, which nulls around the epoch of first light as the HI regions transit from absorption to emission. 
Signal along vertical slices is set by the correlation length of the network, a fraction of the Hubble scale. Fig. 2(a) shows an 
idealized situati on where s tructure is seeded only by strings, in the absence of foreground noise. Contours of 0.5 and 1.0 are 
shown. In Figs. 2(b) 2(c) Gfi = 4 x 10 -7 , representing a system primarily seeded by inflation, but with a 10% contribution 



from strings. Fig. 2(b) includes no sky noise, and has contours representing 0.001 and 0.002. Fig. 2(c)| includes noise, and 



assumes ten year's observation from a square kilometre array (D = 1km, x — 10); recall that x = t b E f covcr D . Contours of 
5 x 10~ 6 and 2.5 x 10" 6 are shown. These plots primarily serve to indicate where signal resides in the observing volume. 
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plot the cross correlation signal-to-noise distribution for 
this case in Fig. 2(a) When only strings are present and 



noise is disregarded, as we do here, the overall normal- 
ization, the string tension G^t, drops out of the signal-to- 
noise calculation. In this imaginary Universe, the signal 
is so strong that an observation at each redshift would 
generate a detection ((S/N)(z,£,Az) > 1). 

In the real Universe, the inflationary perturbations 
overwhelm those from strings, as we see in Fig. |2(b)] even 
without including the detector noise. No single measure- 
ment now has significant signal-to-noise ((S/N) > 1). 
Instead, we must sum over many independent redshift 
bins in order to accumulate a significant signal, giving 
a statistical detection. The case presented in this figure 
represents the best-case scenario: no noise and a max- 
imal, ten percent contribution to primordial anisotropy 
from strings. The signal is similar in shape to, though 
lower in magnitude than, the idealized case of a system 
seeded only by strings. The signal is suppressed more 
significantly on small angular scales, since the inflation- 
ary auto-correlations monotonically increases for £, while 
the string-induced signal peaks at £ ~ 60. Although any 
particular observation in the z — £ plane will result in a 
contribution (S/N) 2 ~ 0(1O~ 3 ), the statistical benefit of 
observing roughly 50 multipoles (2 < £ < 50, weighted 
by a factor of 2£+ 1) in roughly 50 redshift bins (observa- 
tions with a 1MHz bandwidth in the range 20 < z < 120) 
gives rise to a total (S/N) ~ ^/(50) 3 10- 3 ~ O(10). Sum- 
ming our simulated data over all £ and over 70 < z < 20 
reveals a (S/N) ~ 3. However, this does not include ex- 
perimental noise, and corresponds to a survey-strategy 
parameter of x = oo . 



boost our signal significantly once the physics of reion- 
ization are better understood. 
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Fig. 2(c) shows the remaining signal once the sky tem- 
perature is incorporated. These data were generated as- 
suming an all-sky survey from a single dish of diameter 
lkm observing for ten years (x = 10); note that sky tem- 
perature completely overwhelms the signal for redshift 
z > 50. The estimated signal-to-noise remaining in our 
analysis volume is (S/N) ~ 0.28. 



V. OBSERVATIONAL OUTLOOK 

Finally, let us look at Fig. |3j which summarizes the de- 
tectability of our signal for various string tensions, G/i, 
as a function of our observational parameter, x. It is 
clear that the noise induced by the large sky tempera- 
tures severely limits the possibility of observing a string 
network through cross-correlation. However, the calcula- 
tion reported here does not exhaust the possible locations 
of cross-correlation signal from strings. For instance, our 
estimate is restricted to the domain z > 20 that we are 
certain will not be affected by reionization. Our present 
ignorance about the physics of reionization means that 
we cannot trust our calculation into the z < 20 range. 
This is a major impediment, though: lower z is the red- 
shift range least affected by noise, since noise goes as an 
inverse power law with frequency. Hence, we can hope to 
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FIG. 3: Signal-to-noise as a function of observing strategy 
(represented by x = t b B j 'cover D 2 , in units of years x km 2 ) 
and the string tension G/i. Contours of S/N = 1, 2, and 3 are 
shown. 

We have also only looked at the CMB-21 cm cross- 
correlations in this study. Although they are likely very 
small, we are currently analyzing contributions from the 
0(1O 3 ) 21 cm-21 cm cross-correlations. These will also 
give some boost to our signal. 



VI. DISCUSSION 

By incorporating the segment model of strings into the 
CAMB_sources Einstein-Boltzmann code, we have calcu- 
lated the cross correlation between the 21 cm brightness 
temperature and the CMB. Our chief result is plotted in 
Fig. [3j which represents the signal-to-noise available in 
CMB-21 cm correlation studies given an all-sky survey, a 
string tension G/i, and our observing strategy, which we 
characterize by a single integration time and telescope 
properties parameter, x. Unfortunately, noise from the 
sky temperature overwhelms the signal we are calculating 
for observationally-allowed cosmic string tensions. 

However, there is room for more work. We have thus 
far only included the signal from the CMB-21 cm corre- 
lations for z > 20, where we do not need to know the 
physics of reionization. Since the signal should extend 
until the epoch of reionization z ~ 10, we may eventu- 
ally be able to include O(10 2 ) more redshift bins and 
O(10 3 ) data points in the z — I plane, once we can accu- 
rately model reionization. Another contribution to our 
signal that was neglected in this study can be found from 
a cross-correlation study between different redshift bins, 
and this signal will be addressed in future work [5T] . It is 
possible that we will find a detectable signal when these 
extra sources of data are included in the analysis. 

Though our calculations have not predicted a de- 
tectable signal, we reiterate that any detection of a cross 
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correlation between the 21 cm radiation from z > 20 
and the CMB would be a clear sign of new physics. 
Active sources, such as cosmic strings, are a promising 
candidate for the sort of new physics that could gen- 
erate this correlation. However, the fact that inflation 
very likely provides the dominant source of perturba- 
tions in the Universe makes detectability of any signal 
difficult: the inflationary perturbations must necessarily 
enter as "noise" for any study like this, and are typically 
much larger than any alternate sources of perturbation, 
at least prior to recombination where measurements of 
the CMB limit any non-inflationary perturbations. One 
way around this would be to contrive some active source 
that was suppressed prior to recombination which subse- 
quently sourced a larger fraction of perturbations in the 
dark ages. On the other hand, a non-detection of such a 
cross correlation would constitute a strong constraint on 
active sources with such unusual properties. 
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Appendix A: 21 cm vector modes 

In the absence of active sources vector perturbations 
quickly decay and, for that reason, are frequently ignored. 



In particular, they are not included in the 21 cm exten- 
sion to CAMB |20j . Cosmic strings, however, actively 
source vector modes and provide a significant contribu- 
tion to CMB and 21 cm anisotropics. A large effort 
was spent incorporating their effects into our numerical 
code [21] , and we elaborate on some of the details here. 

Applying the harmonic decomposition described in |23j 
and Eq. ^ to (fluctuations in) the photon distribution 
function / provided in |20j gives scalar m — 0, vector m — 
±1 and tensor m — ±2 contributions to the expanded 
field: 



Sf(r},x,e,h) = 




The photon fluid is characterized at some conformal time 
r\ and position x as photons coming from some direction 
n with an energy e = a{rf)E2i- The quantity e repre- 
sents the Doppler shifted energy of a 21 cm photon, E21, 
stretched by the scale factor a(i]). k is the Fourier con- 
jugate to the spatial coordinate x, we denote k = \k\. 

Expressions for the harmonic multipoles F> m ^ are for- 
mally obtained by integrating the Boltzmann equation 
for 21 cm photons along the line of sight [18]. The scalar 

contribution has been reported in |20] , while for vec- 
tor modes we find 



F^(vA,e, k) 



-e- T - (/(e) • [«« + 

— i d 



dr\f c e Tc /,ine(e)v 



(i) ji(kx) 



(A2) 



V(^ 1! - « (1) ) + -v {1) + (r T - 1> (1) 



) e k Xe 



-e- T </(e) 



T.-T, 



7 !' = 



l'=2 



Jn ' v : " i v \, kxT 



r 



for I > 1, where a prime denotes a derivatives with re- 
spect to the argument and a bar denotes a background 
value. Here, t c is the Thomson scattering optical depth, 
t £ is the optical depth to 21 cm radiation, is the vec- 
tor component to the baryon velocity, is the velocity 
(dipole) of the photon distribution, / is the background 
photon distribution function, T 7 is the photon tempera- 
ture, T s is the spin-flip temperature of the neutral hydro- 



gen and determines how much the hydrogen is emitting, 
r T = r e e~ Te /(l — e~ T ') : H is the conformal Hubble pa- 
rameter, Xe = Va — is the conformal distance along the 
line of sight, and ji(kxe) is the spherical Bessel function. 
Quantities subscripted with e are evaluated at the con- 
formal time rj e . The are the angular moments of the 
Fourier expansion of the CMB temperature anisotropy. 
All quantities are discussed at length in [2"U] . 
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The 21 cm power spectra referenced in the main 
text were given in terms of the brightness temperature 
T5. This relates to the photon distribution function / 
through the relation 



T h = 



Eobshpf 

2k b 



(A3) 



sion, hp is the Plank constant, and kf, is the Boltzmann 
constant. Thus, anisotropics in the 21 cm brightness tem- 
perature emitted at a rcdshift z are given as 



8T*{3, h) 



Eobshp 
2k b 



(A4) 



where i? b s is the observed frequency of the 21 cm emis- 
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